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I. INTRODUCTION 


The worldline formalism [IH^ is a field theory technique with various computational advantages such as the reduc¬ 
tion of the number of diagrams in perturbative expansions. It is particularly powerful for amplitude computations in 
external backgrounds as string-inspired techniques become highly efficient in this case. Moreover, fluctuation- 

induced quantities such as quantum actions, energies, and forces can also be efficiently computed by means of Monte 
Carlo simulations of the worldline path integral. Even fully nonperturbative problems can be tackled within the 
worldline formalism [mi. 

The numerical method has given unprecedented access to background problems with nontrivial spacetime depen¬ 
dencies in QED QCD [20,[2l[ fermionic systems [11, as well as the physics of Casimir forces l23-[30l : analogous 

analytical results for QED have also come within reach using worldline instanton approximations [sj® Important 
advantages of the Monte-Carlo method are that spacetime is not discretized but remains continuous in the algorithm, 
and renormalization can be performed on the level of finite and numerically controllable quantities. In addition the 
worldline formalism offers an intuitive picture of quantum fluctuations, as the worldlines can be viewed as spacetime 
trajectories of the fluctuating particles. The numerical method has also been applied to Minkowski-valued 2-point 
functions for external legs [^ . 

In the present work, we demonstrate that the technique can also be extended to composite operators, using the 
induced energy-momentum tensor (EMT) as an example. The quest for such a method was first initiated in [s^. The 
new difficulty arises from the fact that standard tools such as point-splitting lead to expressions involving the 2 -point 
function of the fluctuating field (internal line). As shown below, this can nevertheless be dealt with in a comparatively 
straightforward fashion by using open in addition to closed worldlines. 

The present conceptual and in large parts technical work is motivated by the computation of energy-momentum 
tensors for Casimir configurations. This problem is paradigmatic as the EMT sources the gravitational field. Since 
induced Casimir energies are typically negative (if associated with an attractive force) [d^-j^ , information about the 
EMTs for Casimir configurations is a testbed for conceptual considerations of the interplay of quantum field theory 
and gravity, see, e.g., (45l - l47l | for a discussion of how the Casimir energy “falls”. 

Generically, many fundamental investigations of dynamical spacetime properties start with assumptions imposed 
on the properties of the EMT on the right-hand side of Einstein’s equations. Most prominently, in order to exclude 
exotic phenomena [dsl - f^ . energy conditions (ECs) on the properties of the EMT are imposed. For instance, certain 
“mechanisms” for superluminal travel can be excluded if an energy condition of the form > 0 is imposed, 

where is the tangent vector of a geodesic 7 . If is a null vector this condition is called the null energy condition 
(NEC). The Casimir effect violates some energy conditions such as the NEC [H^ - f^ . It can therefore not be used 
to rule out superluminal travel. In fact, superluminal (but still causality preserving) phase velocities are known to 
occur in the parallel-plate Casimir configuration [ssl - l^ . The Casimir configurations studied so far at least obey the 
weaker averaged NEC (ANEC). For a further discussion of this topic, see [^. Iba - I^ . In the present work, we use the 
violation of the NEC by the Casimir effect as an illustration for our new computational method. For conventional 
computational strategies to determine the induced EMT, for example, by means of mode summation or expansion, 
image charge methods, or similar techniques, see, e.g., [69l - [7^ . 

This paper is intended to be a manual for the numerical computation of composite operators using the worldline 
formalism. It is organized as follows: in Sec.|TT]we apply the worldline formalism to composite operators, specifically 
the EMT of a real scalar field obeying Dirichlet boundary conditions (BCs). We discuss the differences between 
the formalism for composite and non-composite operators or functionals like the effective action. In Sec. m we test 
our numerical worldline algorithm by calculating components of the energy-momentum tensor for a single plate with 
Dirichlet BCs analytically and numerically. The numerical calculation of the EMT and of the NEC is presented for the 
parallel plate configuration in Sec. IIVI We compare our numerical results with known analytical results and discuss 
the arising systematic and statistical errors. We conclude with a summary of our results and an outlook on future 
worldline calculations. 


II. WORLDLINE FORMALISM FOR THE ENERGY-MOMENTUM TENSOR 

Composite, or local, operators in quantum field theory are local products of field operators and their derivatives. 
Being distribution-valued, their product at the same spacetime point may be ill-defined and give rise to divergences. 
Such operators can nevertheless be used for calculations after regularizing the divergences, for example, by point¬ 
splitting, function, or dimensional regularization. In the following, we compute the vacuum expectation value of the 
energy-momentum tensor operator of a scalar field; a composite operator constructed from the scalar field operator 
and its derivatives. 
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Our presentation of the basic formalism follows the one of with only a few differences. We study a quantum 
scalar field ^{x,t) that is a map from a domain V oi d + 1-dimensional Minkowski spacetime onto the 

linear space of self-adjoint operators on the Fock space J- 

l>(x,t): ^ V ^ T. (1) 


t) has mass m and is minimally coupled to a static classical background potential a{x). Later on a{x) will 
impose boundary conditions on the fluctuations of $(x, t). We use the d + 1-dimensional Minkowski metric g^i, with 
the ’’mostly minus” signature (-I-, —). From the Lagrangian density operator 

C = ^df,^{x,t)d>"^{x,t) - i {rr? + a{x)) l>(x,t), (2) 

we derive the equation of motion for $(a:, t), 

{dxd^ + + a{x)) ^{x,t) = 0, (3) 

as well as its canonical energy-momentum tensor operator ($ := $(x, t)) 


Tfj,^{x,t) = 


dC 






= 9^$9i,$ - (^9^$9 a$ - {w? + a 




(4) 


The EMT operator in Eq. O contains products of field operators at the same spacetime point that can lead to 
divergences. We regularize them by a point-splitting procedure in the spatial components x, i.e., 


t) —>■ $(£, t) 

The point-split EMT operator is then ($' := 

T^^{x,t) = Jim 


9a$(x, t)9/3$(x, t) ->■ 9a$(x, t)9^$(x', t). 


- -g^, {d^di -m^- a{x)) 


(5) 


( 6 ) 


Our final goal is to compute the effects of boundary conditions on the energy-momentum tensor. For the sake of 
convenience, we restrict ourselves to computing those EMT components that contain information relevant for the null 
energy condition (NEC) in the z direction, where the z coordinate is always the dth coordinate of our spatial vectors 
X = (a;^,..., x‘^) = {x^, ..., z). This null energy condition is given by the vacuum expectation value of the projection 
of the EMT onto a null curve with the null tangent vector = (1,0,..., 0,1), 


0 < = (foo{x,t)+t,{x,t)) = Tooix,t)+T,Mt), (7) 

where (• ■ ■) denotes the vacuum expectation value. All our calculations generalize straightforwardly to other EMT 
components. In order to compute the vacuum expectation values of Eq. m, we expand $(T, t) in terms of spatial 
eigenmodes and corresponding creation and annihilation operators and dj, 

where a discrete notation is used for both discrete and continuous parts of the spectrum for simplicity. The 'ipj (x) are 
dehned as normalized eigenmodes of the Laplacian in the presence of the potential cr(x) 

+m^ + a(x)^ ^j(x) = E^'tPjix). (9) 

It is convenient to parametrize the energy eigenvalues also in terms of momenta kj according to Ej = kj+m^. Then, 

the eigenvalue equation reads (—+ a{x))tpj{x) = kjtpj{x), which is a Helmholtz-type equation. The corresponding 
Green’s function equation. 


(x) dj + V’y (^) e 


-iEit 


( 8 ) 


-I- a{x) — k^^ Ga{x, x', k) = S (x — x '), 


( 10 ) 
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is solved by the spectral representation 


G^{x,x',k) = 


il)j{x)il}*{x') 
k'i — k'^ — ie 


( 11 ) 


We aim at expressing Eq. 0 in terms of Ga{x,x', k). We therefore insert unity in the form 

1 \ 


CJCJ tXJ / 

1 = J dk 2k S [kj — k^) = J dk 2k lim^ — Im | 

n n \ 


0 0 
and exchange the k integration and the j summation. This yields 


kj — k'^ — ie J ' 


Too{S,t)l = ijrnjdk ^ • (V + V')] 

0 


2Ek 

A: / 1 „ , 1 


^j(x)'ip*{x') 

/j 2 _ ^2 _ 

3 ^ 


T„(f, t)l = faj ik - 1- _V . (v + V') j Im^ 


'lpj{x)'tp*{x') 


k'^ — k'^ — ie' 

3 J 


( 12 ) 


(13) 


The sum inside the imaginary part is just Gcr{x,x' ,k). In general, instead of unit y a decaying exponential 
must be inserted in order to construct local counterterms for renormalization (cf. [^). A then serves as a cutoff 
for large momenta k. This is, however, not necessary in our calculations because we are going to evaluate the EMT 
only at points for which ^{x) = 0, so that all local counterterms, that is, all counterterms that depend on a{x), are 
automatically zero. Furthermore, since the term il)j(x)il)*{x') is in general complex, pulling it into the argument of 
the imaginary part generates an additional term that is proportional to lva'ipj{x)'ip*{x'). This term vanishes in the 
limit X ^ x' and hence will be ignored. The effects of boundary conditions imposed on the field fluctuations by (j{x) 
are described by the difference between the EMT with non-vanishing potential and the EMT with a{x) = 0. So far, 
we have left the potential (j(x) arbitrary to emphasize that our calculations are independent of its specific properties. 
We can therefore repeat all the above steps with a vanishing potential. The only changes that occur are the mode 
functions in E*4- ®. The corresponding Green’s function is Go{x,x',k), defined by Eq. (fTUl) for a{x) = 0, 

rather than Ga{x,x', k). We switch to the common point variables 


^ X + x' 
'■= 2 ’ 



(14) 


and with G{x, x',k) = Ga{x, x',k) — Go{x, x',k), evaluated a.t x = Xcp + A and x' = Xcp — A we have 


^Oo(^cp;t) — 


lim 

A ->-0 




ImG(a;, x', k), 


T„(x„t) = ItaJdG (ji- (8f„ -81) 

0 



lmG{x,x', k). 


(15) 


Equation (1151) is independent of any specific mode expansion of $(x, t). It only depends on the Green’s function 
G{x,x',k), which is representation-independent by definition. Therefore, any method for computing G{x,x',k) can 
be used at this stage, e.g., an optical approach in [t^- The subtraction of Go{x,x',k) also removes the divergent 
terms that are independent of the potential (t(x) and normalizes the EMT such that it vanishes for (t(x) = 0. 


A. The worldline representation of G{x,x',k) 

For simplicity, we consider here only the massless case m = 0. Our calculations can straightforwardly be carried 
over to the massive case, for details, see App. In order to express the Green’s function in the worldline formalism, 
we interpret G(x, x', k) as the matrix element of an operator G{k) and Eq. (j9|) as a quantum mechanical Schrbdinger 
problem whose Hamiltonian is E[ = —-I- cr(x). Then G(x, x', k) corresponds to a quantum mechanical propagator, 
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Fourier transformed to energy space, from which the free motion has been subtracted. Hence, it can be written in 
position space in terms of a propertime representation: 

OO 

G(f, x', k) = (x'l g{k) \x)=i J ds (x'l |x). (16) 

0 

The matrix element in Eq. (1161) is the quantum mechanical transition amplitude of a fictitious particle moving from 
X at the fictitious time r = 0 to x' at r = s with a Hamiltonian H = — + cr(x). The corresponding Feynman path 

integral in position space is then straightforward to find. Next, we perform formal Wick rotations in both the s and 
k planes, such that s = —iT and k = —ikE, which are consistent with causality. This casts G(x, x', k) into its doubly 
Wick-rotated form 


OO x'=S{T) T _ 

dT e“^^® / 'Dx{t) e 0 * 

0 :e—;e( 0 ) 

The variables s and T are called Minkowskian and Euclidean propertime, respectively. Both describe the time 
evolution of the fictitious Schrodinger problem, but neither is a physical, measurable time. An analogous propertime 
representation for the effective action has been used in previous calculations of effective interaction energies for the 
Casimir effect and similar boundary configurations [ll|, [2^, [2^ . The worldline representation of the effective action 
contains, however, a path integral over closed loops, whereas for G(x, x',fc), open worldlines running from x to x' 
must be computed. The implicitly normalized path integral in Eq. dlB gives, for the free path integral, the standard 
free propagator. 


- /dr (t(5(t)) 

e 0 -1 


(17) 


x'=x{T) 


Vx{T)e 0 


rr=rr(0) 


(47rT) 


d/2 


with A = 


(18) 


From this normalization one derives the shorthand notation for the path integral expectation value, the worldline 
average of an arbitrary operator O [x(r)]: 


{0[S{t)])x,x' 


x'=x{T) T 

/ Vx{T)e 0 G [x(t)] 

x—x{Q) 


g-A^T 
(4^y)d/2 ■ 


(19) 


Pulling out the imaginary part prescription in front of the k integral in Eq. m. we can perform both Wick rotations 
for the components of the EMT. Details can be found in App. |Al We finally arrive at the manifestly real worldline 
expressions for the normalized vacuum expectation values of Tqq^x,/) and Tzz{x,t), 


^00 {Xcp ; l) 


^ZZ (Xcp ; 


lim 


dT 


/A \ — I A m 

(47r) 2 A-j-oJ 4T 2 
' ' 0 



T 

-/dTcr(x(r)) 

e 0 -1 


X cp, ^ 



-/ dTG{x{T)) 

e 0 -1 


□ ,A 


( 20 ) 


The expectation value of the path integral in Eq. dini) is defined by Eq. (HU) with the operator O [x(r)] = 
exp /J"dTcr(x(T))^ — 1. It is an average over worldlines x(t) that start at x = x(0) and end at x' = x(T), 
and that are weighted with a Gaufiian velocity distribution. It is convenient to switch to coordinates such that for 
X —>■ x' all paths are fixed with respect to one common point Xcp = {x + x') /2 (cf. Fig. [T]). We therefore call these 
worldlines common point loops or lines. The remaining path is written in terms of a dimensionless unit worldline y(t) 
with t G [0,1] and t = T t. The unit line y{t) itself can be written as the sum of the classical path from x to x' and 
a deviation Z{t) that obeys 2^(0) = 0 = 2(1). We choose the origin of the y coordinate system to be the point Xcp, 


x ( t ) = x{Tt) = Xcp + VTy{t) = Xcp + Vt (^(1 - 2t)5 -\- Z{t)^ . 


( 21 ) 
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FIG. 1. Schematic depiction of the rescaling from the worldlines x(t) to the paths x^p + y/Ty(t) with t = Tt and the unit 
worldlines y{t). 


In the limit 5 —)■ 0 the worldline y{t) closes, it becomes a loop y{t) —>■ Z{t). We have set above 


A x — x' 

Vt ~ 2Vt 


( 22 ) 


The VT dependence in this relation is crucial for the computation of the Aj;-derivatives in Eq. (Uni- Rescaling the 
worldlines, x{t) —^ y{t), makes the weight factor of the resulting y path integral independent of T. This allows us 
to compute the expectation value by generating only one ensemble of unit worldlines y(t), which are all defined with 
respect to the same coordinate system. The worldline average is now calculated by replacing the path integral with 
a sum over a finite number N of unit worldlines that are themselves approximated by a finite number Uppi of points 
per line yi with i G {0,... ,nppi}. The points y) are random numbers distributed according to the Gaufiian weight 
factor exp(— dty^/4), being a discretized realization of an open line with the vector 2S connecting its endpoints. 
As a result, the expectation value of an operator O is written as an average, 


1 ^' 

^ ~ ^ * G {!) • ■ •) fippz}- (23) 

Xap, V 

The common point paths yi are conveniently generated with the d loop algorithm [llj |. which also works for open 
worldlines. 


B. Dirichlet constraints on dT> 

The entire formalism that we have outlined so far works, of course, for arbitrary static background field configura¬ 
tions u(x). In our calculations, we use 


a{x) = X / dll5{x — xdv)i Xdv & dV, (24) 

dv 

where dS is a surface element on dV. We recover Dirichlet BCs in the limit A —>■ oo (cf. [zllii)- The subtraction of 
the vacuum Green’s function in Eq. (HID already removed all divergences independent of a{x). Therefore, G{x,x',k) 
can only diverge at points for which a{x) ^ 0 and the energy-momentum tensor in Eq. (1201) is finite on (D, where tr 
vanishes. All remaining divergences are located on the boundary dT) and can be related to the infinite amount of 
energy necessary to constrain 4>(al, t) on all momentum scales to fulfill the Dirichlet boundary condition on dT). 
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The parametrization Eq. (1241) and the subsequent Dirichlet limit greatly simplify the worldline average because now 
we have 


1 



-1 


3- cp 7 ^ 


-1 

0 

0 


if Xcp + intersects dV 

otherwise 


S {xcp + VTy{t)) 



(25) 


Equation (1^51) states that only paths x^p + VTy{t) which violate the boundary conditions lead to deviations from the 
trivial vacuum and thus contribute to the expectation value. We call the function S the intersection condition. It 
gives a geometric description of how and for which values of T the worldlines intersect the boundary. As long as the 
start and end point of a given worldline are not on the boundary (alcpi A ^ 922), the intersection condition will always 
determine a minimal non-zero value for the propertime T^m > 0 for which the worldline first intersects 922. Tmin 
denotes the minimum propertime that is necessary for a worldline to propagate from the start to the end point and 
intersect a boundary in between. For any T < Tmin deviations from the straight line between x and x' are typically 
strongly suppressed. Only for sufficiently large propertimes T > Tmin does the diffusive Brownian motion process, 
described by the path integral, create sufficiently large random detours that can intersect 922. In the propertime 
integral, Tmin serves as a lower bound and removes the divergence for T —>■ 0. From a physics point of view, Tmin acts 
as an ultraviolet cutoff because small propertimes correspond to large momenta. The intersection condition S may 
also provide an upper bound Tmax- This value Tmax would then mark the propertime for which the worldlines no 
longer intersect the boundary. This needs to be considered especially for configurations where the boundary consists 
of compact objects (see e.g. [^1. 


C. Compact expressions of Too{xcp,t) and Tzz{xcp,i) for worldline numerics 


The EMT components are completely finite on 22\922 by virtue of Eq. (l25l) . We may therefore decompose Too(xcp, t) 
and Tzz{xcp,t) further and study the resulting, more compact, terms one by one. Toward this end, we parametrize 
the intersection condition as S {xcp + VTy{t)) — y/TM. — 1, where the function A4 has inverse length dimension 
and describes the geometrical conditions for a worldline to intersect the boundary 922. It depends on Xcp and <5, and 
through the latter also on T. This T dependence is not explicitly known but it vanishes as A —>■ 0. In general, A4 is not 
a smooth function of these parameters. In fact, it can be non-differentiable in either variable. Since A4 depends on the 
shape of the boundary 922 for the setup under consideration, its functional properties must be investigated each time 
anew. In addition, there can be several ways in which the intersection condition S may be written. Whereas different 
parametrizations are equivalent, they can exhibit very different properties during numerical evaluation. Thus, the 
following calculations should be understood as primarily formal. First, we assume the function Ai to be differentiable 
twice in both Xcp and <5. We additionally assume that all limits can be evaluated, and that M. only determines a 
lower bound Tmin on the propertime integral. In this case, we have Tmin = 

We now insert Eq. (|2^ into Eqs. (|2Q|) and find 


^00 (^cp; 0 


Tzz 5 0 


OO 

1 r dr 

47 (4^T)^ 

0 ' ^ 

OO 

1 r dr 


(c - - V%) e-^ {q[Vtm - l) 


(26) 


In order to evaluate the propertime integral, we need to interchange several limits, which requires justification. First, 
we note that all expressions we deal with are finite by construction due to the Dirichlet constraint Eq. (EH). Since 
this is valid for all A, the limit A —^ 0 may be interchanged with all other limits except 9 a^. Furthermore, since we 
approximate the worldline average by a finite sum, this average can be exchanged with other limits and be conveniently 
computed last. We may also formally interchange propertime integration and differentiations because T, Xcp and A 
are independent variables. There is, however, the intersection condition S = y/TAi — I, which depends on all three 
variables. The function Ai also depends on the combination <5 = K/y/T. As a consequence, the derivatives of 5, or 
more specifically of AI, must be computed before the propertime integration (cf. Sec. |B]). Picturing the worldlines 
as paths in space helps examine the situation: we need not only compute the intersection condition itself but also 
its derivatives, that is, we need to determine how the intersection is altered if Xcp or A are changed. A derivative 
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with respect to Xcp can be viewed geometrically as moving the complete worldline through space without changing 
its shape. By contrast, the A derivative corresponds to opening and closing the worldline at a fixed point Xcp in 
space, changing its shape in the process. The required order in which these manipulations of the worldline expressions 
should be performed is then: 

f. compute the derivatives of 0 (^VTJ\4 — with respect to Zcp and Sz = I^z/VT, 


2. let A —0, that is, 5 —>■ 0, 

3. perform the propertime integration, and 

4. average the expression over all worldlines in the ensemble. 

With these consideration at hand, we write Too(xcp,t) = 2oo{a?cp,0|/ + 2bo(^cp,i)|jj where we define 

1 1 


^Oo(^cp ;0 r ■ — 


(47r)“2“ d + 1 




^CV - 0 


Toq{xcp, t) \ 


1 1 


1 


_y72 

2(4^)4id-l\ "P /x„„A=o 


With Tzz{xcp,t) we proceed accordingly and write it as a sum of four terms: 


TzAxcp,t)\^^ 2 d-1 


2 (47r) 2 d 1 

Tzzix,p,t)\j^ := - 


d^Zz{_^Cp-! T T ■- 


1 1 




^ (47r) id — 1 \ ^ / Xcp,A—0 


Xcp,A — 0 ’ 

Xep,A=0 ^ -Too{Xcp,t)\j , 
— Too{Xcp:t') 


II ■ 


(27a) 

(27b) 

(28a) 

(28b) 

(28c) 

(28d) 


The fourth term, Tzz{xcp,t)\j^, comes from acting with d\^ on the exponential e ^ Another term, the mixed 
term (9^^ e“^^/^)9A^ (■ ■ •) vanishes as A —>• 0. 

We note that, with the help of the above decomposition into compact terms, the null energy condition along the z 
axis Eq. © now reduces to computing 

0 < roo(x,t) +Tzz{x,t) = Tzz(xcp,t)\j^ +Tzz{xcp,t)\j^^. (29) 


III. THE ENERGY-MOMENTUM TENSOR FOR A SINGLE PLATE 

The numerical computation of Tqq (a;, t) and Tzz (a:, t) in d = 2 and d = 3 space dimensions in the case that dT> is 
a single d— 1-dimensional surface, i.e., a plate, is our first proof-of-principle example. The plate imposes Dirichlet 
BCs on the fluctuations of $(x, t). It is placed at z = 0 such that its normal is the z axis. The single Dirichlet plate 
configuration is also sometimes referred to as the perfect mirror. 


A. Analytic calculation for a single plate 

Before we use worldline numerics, we compute the EMT for the single Dirichlet plate analytically. Eor that, we 
use Eq. (ITSl) and compute G{x,x',k) by solving the equation of motion Eq. © for different boundary conditions. 
Denoting the BCs by a superscript cr, we must solve 

(_V 2 -p 2 ^^;(f) = 0 . 


(30) 
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Toward that end, we decompose any x G 2? in the d — 1-dimensional vector x\\ parallel to dT) and the z component 
of X. The solution of Eq. O in the half space z > 0 with Dirichlet boundary conditions at z = 0 is then 


V'p (®) = ^ (exp {ip- x) - exp (ip • 


(31) 


where x = {x\\,z') and x = (x||,—z). We assumed that there are only outgoing waves at spatial infinity (Som- 
merfeld radiation condition). In the same manner, the free solution without boundary conditions is found to be 

i/'p(x) = exp (ip-x). Both solutions are normalized for z > 0 and z' > 0, / 'ipp{x) ipp* {x') = 6^ {x — x') = 

f '’Pp(x) 'ipp*(x')- The Green’s functions Ga{x,x',k) and Go{x,x',k) are now computed according to the spec¬ 
tral representation Eq. (HU. We perform the momentum integration in polar coordinates and find for G{x,x',k) = 
G„{x,x',k) — Go{x,x',k) (see also [zilzS) 


1 

2tt 

1 

47r 


G{x,x',k)\‘^ ^ = — (lCo{—ik\x — x'\) — /Co(—ifc|x — x 
ZTT V 

_ z\x—x'\ f^ik\x—x I \ 1 f^ik\x — x'\ 

G(x,x',fc)r=" " '" ' 


ZTT 


(32) 


47r lx — x'l 


The /Co are modified Bessel functions of the second kind (Macdonald functions). 

We can now solve Eq. (USD analytically and even use the decomposition of Too(xcp,t) and Tzz{xcp,t) that we 
developed in Sec. Ill Cl Denoting the distance from the plate with Zcp, we find for the EMT in d = 2 


^ 00 (^cp; 0 |/ — 


1 1 
327r z^3 ’ 
1 1 


IGtt z 


3 ’ 
cp 


Too{Xcp,t)\jP — 267rz3’ (33a) 

T'2:z(xcp,/) I = 0. (33b) 


And in the case of 3 spatial dimensions the values for Too(xcp)/) and Tzz{xcp,t) are 


'^Ooi^cp-) 


1 1 

^Oo(^Cp5 0 

1 d—3 

3 1 

/ 

327r2 z4 ’ 

// 

327r2 z4 

^zzi^cpi i) 

\d^3 
\ la 

3 1 

327r2 z^ ’ 

^ZZ (^cpj 1 

d^3 

0. 


(34a) 

(34b) 


We note that Tzz{xcp,t)\p^ and (/^^(xcp,/)|cannot be calculated separately in a direct manner since we used the 
functional structure of the worldline representation of G{x,x',k) to define these functions. Despite that, we can 
always compute Tzz{xcp,t )\from the sum Tzz{xcp,t )\using the fact that Tzz{xcp,t)\p^ = —Too(xcp,/)|^. We 
thus have 

Tzz{Xcp,t)\j^ = Tzz{Xcp,t)\jf^ = 

According to Eq. (l29l) the NEC along the z axis is then violated. 


-/bo(^cp51) T Tzz(xcp, /) 


— 1/ (327rzj^p) for d = 2 

— 1/ (IGtt^ Zgp) for d = 3 


These are the same values for the NEC which were derived in 


(36) 


B. Worldline calculation for one Dirichlet plate 


The first step in all our worldline calculations is determining the intersection condition S = y/TM — 1. For the 
single plate setup, this is easily done from Fig. [2l The worldlines y{t) start at the point Zcp and intersect the plate at 
z = 0 for all T that fulfill y/Ty- + Zcp < 0. We call the z component of the point on the loop that is closest to the 
plate p_. It is negative for our choice of coordinates in the setup Fig. [5] Thus we find for the 0 function 


y/Ty^ + Zcp <0 


0 


(37) 
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This corresponds to a minimal propertime Tmin = ■, where (—?/-) is the positive distance that measures the 

extension of the loop towards the plate. Tmin only depends on the z component of yit) because the plate constrains 
the propagation of j7(t) only in the z direction. For this reason, and because the worldline distributions factorize 
with respect to their position space components, we only need to calculate 1-dimensional loops. Furthermore, only 
one point of every loop, the point ?/_, needs to be found. In this case, we have Tmax oo, he., if the intersection 
condition is fulfilled for Tmin-, it will be fulfilled for all T > Tmin- 

The 1-dimensional open line y(t) depends on (5 and thus on the ratio AzIVT- Being an extremal point on this loop, 
y- := y{t-) carries the same dependence. The defining equation for the extremum is (cf. Eq. ([Si]) ! 


y- ■= y{t = t-,Sz) 


dt 


dtZ{t)l_-2Sz = 0. 


(38) 


This is an implicit equation for t- which is itself a function of 5z- Consequently, the minimal point depends on 5z 
in an explicit and implicit way, that is y_ = y{t-(6z),Sz)- In order to compute ToQ{xcp,t) and Tzz{xcp,t), we insert 
M = {—y-)lzcp into Eq. (I27a|) - (l27bl) and Eq. (|28al) - (l28dl) . In Too{Xcp,t)\j, only a power of M must be computed, 
which can straightforwardly be implemented numerically. There are also no difficulties in calculating t) 

and Tzz{xcp,t)\ as A4 is obviously differentiable in Zcp- 

However, M is strictly speaking not differentiable with respect to 6z because the Brownian bridge y{t) is not 
differentiable. This problem is caused by a naive interchange of derivatives and the path integral. Our strategy is to 
cure such problems by identifying a constructive definition of the derivatives in terms of geometric properties of the 
discretized worldlines. Eor the moment, let us write 


Tzz{Xcp,t)\pf^— 

(47r) 


1 {(-y)d(-y y' + d{-y.r-^i-y.y{-y.y) 


2 


with (•)' = ds^{-)- We compute the derivatives with the help of Eq. m and (1551) and find 




<5^=0 


= -(l-2t_) 


5^—0 


S^=0 


= 2 


t-{Sz + hs) - t_{Sz - hs) 
2hs 




(39) 


(40) 

(41) 


The right-hand sides of these equations are, in fact, well-behaved, such that they can be used as a definition of the 
derivatives of y{t) to be inserted into the path integral. The constraint Eq. (1551) allows for an analytical form of (—y_)'. 
For its computation only the time parameter t-{5z = 0) must be determined. A similar closed form of (—?/_)" does not 
exist because the function t-{5z) is unknown. It can only be calculated numerically by using the difference quotient 
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in Eq. dun with <5^ = 0 and sufficiently small hs- The ds^ derivative terms are not necessary for the computation 
of the NEC in the case of the single plate because we know Tzz{xcp,t)\= —Too(a?cp, t)|^ (in general) and from the 
analytical resnlts = 0 (for the single plate). Hence, we have Tzz{xcp,t)\j^ = Too{xcp,t)\j- and we 

know that Tbo(aicp, t) |^ contains no derivatives. We use this as an additional check for computing derivatives with our 
algorithms. The null energy condition along the z axis can consequently be read off from TQo(xcp,t) + Tzz{xcp,t) = 
^zz {Xcp-, 0 I /fj T ^00 {Xcpi I / ■ 


C. Numerical results for TQo{xcp,t) and Tzz{xcp,t) 

We need to compute six different worldline averages, ((—?/_)p) with p G {1, 2, 3,4} and (dg {—y-Y) with q G {3,4}. 
As with any Monte-Carlo method there are statistical and systematic errors. The statistical error arises from the 
finite number of worldlines N in the ensemble. We use the standard deviation of the observables as a measure for the 
statistical error, which decreases as Since the computational costs of our worldline algorithms scale linearly 

with N, the statistical error is readily controlled. 

The systematic errors in our calculations arise from two different sources. The first arises from the discretization of 
the worldlines themselves and is controlled by the number of points per lines Uppi. The continuum limit is defined by 
Uppi —>■ oo. In Fig. [3] we show this by plotting our worldline results for ((—j/_)^') with p G {1,2,3,4} against Ijuppi. 
For increasing Uppi , the numerical data approach the analytical values from below. We fitted our worldline data with 
a function a + b/uppi, disregarding the error bars, i.e. the statistical errors, as they do not vary much for different 
Uppi. The systematic error is then the difference between the analytical value and the fit in the limit nppi —>■ oo. The 
second systematic error is introduced by the difference quotient used to compute the derivatives of the minimal point 
((—?/_)). For the numerator in Eq. (1411) to be non-vanishing, hs should be at least of the order of the average distance 
between points. Only then, the change of <_ as a function of <5 can become visible; for too small hs, t_ simply does 
not vary such that y” —^ 0 arises as an artifact of the discretization. At the same time, hs must be small compared to 
the extension of the entire loop. The variance of the points of the loop (j/^) = 1/6 [2^ gives a rough estimate of 

the extension of the loop. For the average distance between points, we use (jYi=o iVi+i ~ Vi)^ •\/2/ Uppi. 

Both values are only rough measures, they become satisfactory approximations only for large N and Uppi. Combining 
both requirements, we find 


hs = \ i with / > 0(1) / = I ■ with e < 1. (42) 

y nppi 6 6 V 2 

Since Eq. (HI is just an estimate, we determine optimal values for / directly by computing ((—y_)") for different 
worldline ensembles with N = 25 ■ 10^, Uppi = 2^° .. .2^° and several values hs = f s/2/ Uppi. We plot our results 
against / (cf. Fig. |3][7]) and find large errors for small /, especially for / < 1. These values violate our requirement 
Eq. (1421) . For intermediate / there is a plateau region beyond which the values of ((—j/-)") deviate (non-linearly). 
For these / values a linearization is not a sufficient approximation for the derivative anymore. 

We read off optimal values for / for different Uppi from the plateau region in Fig. |3][71 These values are compiled 
in Tab. ID We estimate the optimal / for values of Uppi for which we have not run a calculation, such that the 
corresponding e are approximately 0.1. Those values are consistent with the constraint in Eq. (1421) . In order to satisfy 
/ ^ C)(l) and e ^ 0, our results suggest that computations of ((—?/_)") be only run with Uppi > 2^^ (cf. Fig. |4] 
and[S]). We deduce the systematic errors of averages of different powers of (—y-) from the fits in Fig. |31 Similarly 
we determine the systematic errors of the averages that contain derivatives of (—y-) by computing these terms for 
different values of Uppi and the corresponding optimal value of /. Fig. [8] shows the results of these calculations. For 
small Uppi our numerical data approach the analytical values from below, but are greater than the analytical values for 
large Uppi. We again read off the systematic error as the difference between analytic values and our fit for Uppi oo . 
This concludes our error discussion for the single Dirichlet plate. With the plots and optimal / values from above, 
we are now able to run our algorithms with minimized errors. 

To compute TQQ(xcp,t) and Tzz(xcp,t) with minimal errors we use an ensemble of Ai = 5 • 10^ worldlines with 
nppi = 2^° points per line and choose f = 12 {e = 0.099). The statistical errors are measured by the standard 
deviation, the systematic errors by the difference between analytical values and fits for Uppi —>■ oo. In Tab. HIl the 
worldline results are summarized for these optimized parameters. The statistical errors are well below 1 %, due to the 
large number N of worldlines in our ensemble, and the systematic errors are smaller than 1.5%. They are larger for 
higher powers of (—y_). The errors of the derivative terms are generally larger, for they not only contain the systematic 
error for the point (—y_), determined by Uppi, but also the error due to the linearization of the derivatives, determined 
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FIG. 3. Numerical values for ((—y_)’') with r € {1, 2, 3,4} for Uppi = ... 2^° with error bars showing their statistical errors. 

The solid curves are fits to our data and the dotted lines are the respective analytical values. Systematic errors are the difference 
between analytic results and our fits for Uppi 0. 


by hs- To summarize, the numerical worldline method is capable of computing the relevant EMT components for the 
present case on a satisfactory percent precision level. 


IV. THE ENERGY-MOMENTUM TENSOR FOR PARALLEL PLATES 


Our second proof-of-principle example is the case of two parallel plates with Dirichlet boundary conditions. This 
configuration corresponds to the original setup studied by Casimir for the electromagnetic field . The EMT in this 
setting has been studied in the literature, see, e.g., 0. A.S m the previous section, w^e compute Tqq^x, and 
in d = 2 and d = 3 dimensions. With these terms, we can see if the null energy condition is violated or fulfilled. 
The plates are placed at z = a/2 and z = —a/2 and the z axis is perpendicular to them. The whole configuration is 
invariant under arbitrary translations parallel to the plates. 
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FIG. 4. Numerical values of ((—y_)"} as a function of the difference quotient parameter / computed with three ensembles with 
N = 25 ■ 10^ and Uppi = 2'^°. The three sets of data points have been slightly shifted horizontally to be distinguishable in the 
plot. The solid black line is only for reference and comparison between different figures. 



FIG. 5. Numerical values of {{—y-)'') as a function of the difference quotient parameter / computed with three ensembles with 
N = 25 ■ 10® and Uppi = 2^^. The three sets of data points have been slightly shifted horizontally to be distinguishable in the 
plot. The solid black line is only for reference and comparison between different figures. 


A. Analytical calculation for two parallel plates 


Too (x, t) for two parallel Dirichlet plates can be calculated analytically by applying the method of images and 
using the results of Sec. IIII Al for the Green’s functions. The two plates at z = ±al2 constitute the boundary and 
decompose the domain T) into three disjoint regions. In the outside regions ^ > a/2 and ^ < —a/2, we find the 
Green’s function for a single plate at z = a/2 and z = —a/2, respectively. The Green’s function between the plates 
is given by the method of images as an infinite series of image charges induced on the plates by a point source. For 
arbitrary x and x' between the plates for which \z\ < a/2 and \z'\ < a/2, we introduce Xq := {x\\,z + 2qa') and 

Xq := [x\\,—z+ (2q + l)a). The Green’s functions G{x,x',k) = Ga{x,x',k) — Go{x,x',k) are then 


id=2 1 


G{x,x', k) 

G{x,x',k)\'-^ = 


^ X! {^o{-ik\xq - £'|) - K.Q{-ik\xq - £'|)^ - —K.o{-ik\xq=o - £'1), 


qe i 

gife|xg-S'| \ 2 Qik\xg=o-S'\ 


96 : 


dTT \Xq=o - X'\' 
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FIG. 6. Numerical values of {{—y-)'') as a function of the difference quotient parameter / computed with three ensembles with 
N = 25 ■ 10^ and Uppi = 2^^. The three sets of data points have been slightly shifted horizontally to be distinguishable in the 
plot. The solid black line is only for reference and comparison between different figures. 



FIG. 7. Numerical values of ((—as a function of the difference quotient parameter / computed with three ensembles with 
N = 25 ■ 10® and Uppi = The three sets of data points have been slightly shifted horizontally to be distinguishable in the 
plot. The solid black line is only for reference and comparison between different figures. 


After inserting these Green’s functions into Eq. (I15L we integrate with respect to k. We simplify the result further 
by assuming that x and x' lie on the z axis, that is, a:|| = x!^y The various EMT components can now be written 

in the form of a series C(s, f{zcp, ^z)) ■= (/(-^cp, A^) + q)~^, which is a Hurwitz ( function for Re (s) > 1 and 

Re (/(zcp, A^)) > 0. While we always have /(zcp, A^) > 0, for d = 2 we find s = 1 in Tzz{xcp,t). We solve this 
problem by first acting with the derivatives and which increases the exponent to s = 3. 

In d = 2 the EMT between the plates, \zcp\ < a/2, is given by the components 


^00 i^cp ; 0 


d=2 

I 


Too{xcp,t)\j^ = 

Tzz(Xcp,t)\^l^_^^^ = 


1 

327ra® 

1 

10710 ® 

1 

10710® 

1 

107ra® 







+ C(^3, 

+ Cf3, 



2C(3). 



(43a) 

(43b) 

(43c) 

(43d) 
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TTjppl 

allowed / 

allowed e 

optimal / 

optimal e 

2io 

1.0. 

..1.5 

0.265., 

. .0.398 

1.00 

0.265 

2“ 





1.00 

0.188 

212 





1.25 

0.166 

213 





1.50 

0.141 

214 

2.0. 

..4.0 

0.133., 

, .0.265 

2.00 

0.133 

215 





2.75 

0.130 

216 





3.50 

0.116 

2i" 

2.5. 

..8.0 

0.059., 

, .0.188 

4.00 

0.094 

218 





6.00 

0.099 

219 





8.00 

0.094 

220 

9.0.. 

.25.0 

0.075., 

, .0.207 

12.00 

0.099 


TABLE I. The optimal values for the algorithmic parameters / and e, respectively, in the range Uppi = .. .2^°. They are 

determined from Fig. |3]17|or estimated such that e « 0.1. 


0.97 

0.95 

0.93 

0.91 

0.89 

0.87 

0.85 

0.83 



‘^ppl 


'IT'ppl 


FIG. 8. Numerical values of the derivative terms in d = 2 and d — 3 for Uppi = 2^° ... 2^°, N = 25 ■ 10® and the corresponding 
/ from Tab. |T] The solid curves are fits to our data and the dotted lines are the respective analytical values. Statistical errors 
are given by the error bars, and the systematic errors correspond to the difference between analytical and fitted data in the 
limit Uppi —^ oo. 



analytic value 

numerical value 

Astat [%] 

Asys[%] 

Acomp [%] 

ii-y-)) 

Ri 0.88623 

0.88502 ± 0.00065 

0.073 

0.339 

0.412 

ii-y-r) 

1.00000 

0.99715 ±0.00141 

0.141 

0.500 

0.641 

(i-y-f) 

1.32934 

1.32337 ± 0.00289 

0.218 

0.527 

0.745 

{(-y-)^) 

2.00000 

1.98702 ± 0.00626 

0.315 

0.550 

0.865 

Hdli-y-f) 

0.88623 

0.89190 ± 0.00452 

0.507 

1.467 

1.974 


1.00000 

1.00103 ±0.00561 

0.560 

1.300 

1.860 


TABLE 11. Numerical results (and their standard deviations) for the single plate in 2 and 3 spatial dimensions with an ensemble 
of 5 • 10® worldlines, 2^° points per line and hs = 12 y'2/n^. The statistical and systematic errors are shown explicitly, as well 
as the complete error Acomp- 
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The corresponding 3-diniensional results are 

|d=3 1 


^00 (^cp? I 

^00 (^cp; I 
Tzzi^^cp-i I 

Tzz{Xcp,t)\j^j^j^ 


1 / 

\d^3 

^Cpi f'J I jj 

\d^3 

'^cp^ojli^ 

|d=3 


327r2a4 

3 

327r2a4 

3 

327r2a4 

3 

327r2a4 


2C(4)-C 4,- + ^^ -C 


C 4,- + ^ +C 4,--^ 


4'‘-5 + - 

\ 2 a 

2C(4). 




(44a) 

(44b) 

(44c) 

(44d) 


As before Tzz{xcp,t )\can be computed using Too(a:cp,t)|^ = —Tzz{xcp,t )\The EMT components diverge as the 
distance from either plate goes to zero. The NEC along the z axis between the plates is also violated in the case of 
two parallel plates, 


- [6C(3) + C (3, i + ) + C (3, i - ^)] / (327ra3) for d = 2 

- [4C(4) + C (4, i + 5^) + C (4, i - ^)] / (iGTT^a^) for d = 3 ’ 


Our results agree with those of [t^. 


(45) 


B. The conformal complement oi T^v(x,t) 

The energy-momentum tensor for the electromagnetic field with perfect conductor boundary conditions does not 
between the plates [^. However, the electromagnetic EMT has zero trace because the electromagnetic field is 
conformally invariant. The energy-momentum tensor derived for a conformally coupled scalar field between two 
parallel plates with Dirichlet boundary conditions does also not diverge near the plates but it still yields the same 
total energy and momentum as for the minimally coupled scalar. The difference between both tensors is the conformal 
complement ATiz^{x,t) := {dp,d,z — g^^dad^^) ^{x,t)^{x,t), where ^ = (d — l)/(4d) [t^. Adding the vacuum 
expectation value of the components of this conformal complement operator to our canonical EMT removes the 
divergent terms of the latter. We determine the conformal complement in the worldline formalism along the lines of 
the calculations performed for Too (x, t) ai^dTzz{x,t) in Sec.HIl We find AToo{x,t) = —4^-Too(T, t)|^^ and ATzz{x,t) = 
—4:^Tzz{x,t)\j.^ + 4^Too(af, t)|As a consequence, the conformal or improved energy-momentum tensor Qp,v{x,t) := 
Tp,v{x,t) + AT^^{x,t) is then finite, even constant, but still violates the NEC: 

0oo(*, t) + Qzz{x, t)=- ^ 2C(3) - 2 ^ 2C(3) for d = 2, 

OZTT UZTT ^^0^ 

0oo(^, t) + Qzz{x, t)=- 2C(4) - 2C(4) for d = 3. 

Our results for (T, t) match with those of M where the general case of a massive scalar field with Dirichlet BCs 
in d spatial dimensions is presented. Since we can always construct 0^,^ (T, t) from its canonical counterpart T^,^ (T, t ), 
we restrict our computations to the latter. 


C. Worldline calculation for two parallel Dirichlet plates 


We start our worldline calculation with identifying the intersection condition from Fig. The worldlines have now 
two possibilities to intersect the plates, which means we need to identify the two extremal points y- and y+ of each 
worldline. More precisely, we need to determine which point intersects one of the plates first. For any given line this 

will depend on the position Zcp on the z axis. The intersection condition A4 in 0 ( y/TM. — 1) is then 


XTy. 

VTy+ 


7 < —- 

Zcp \ 2 

Zcp ^ 2 


with Ai = max(A4+,Af-) := max 


{y+) i-y-) 


L 2 


- 2 c 


+ Zc 


(47) 


The condition in Eq. (|47l) is 1-dimensional, so that 1-dimensional worldlines suffice for the computation. However, 
the calculation, especially of the derivatives of A4, is more involved because M is the maximum function AA = 
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{M+ + M- + \M+ — M-\) /2. Due to the modulus, - if taken at face value - is not differentiable at M+ = A4-, 
neither in Zcp nor in S^. As before, we still treat M as if it were differentiable in order to identify a constructive 
definition of the derivative inside the path integral in terms of geometric properties of the worldline. Worldline averages 
of powers of M are straightforward to compute. Averages of first derivatives of M also pose no difficulties because 
dx\x\ = sgn(a::). However, worldline averages of d^M. involve the derivative of the sgn function. Formally, one has 
dx sgn (x) oc 5 (x) whose numerical evaluation under the worldline average is not straightfoward. So far, we have not 
been able to find a suitably smeared-out 5 function adapted to the numerical discretization with controllable errors. 
Instead, we simply approximate 9a;Sgn(x) with a difference quotient. In the discretized numerical framework, there 
is no unique variable choice. For example, one can consider sgn \M.+ — M.-] as a function of one variable A1+ — M.- 
or as a function of the two variables M+ and A4-. The difference quotient of either option suffers from problems. 
In the first version — Ai- can go to zero and will naively generate diverging terms. The second version, while 
avoiding divergences and yielding better results, still has a large systematic error. A more successful strategy arises 
from a third option to view the sgn function as a function of M.+/A4-. This is possible because both and Ai- 
are positive. The resulting difference quotient is thus written as 


9 sgn [M + 


M-] = 9 sgn 



= ^^4+ sgn 


M+ 

M- 


d 


M+ 

M- 


sgn 


^(1 + 0-1 


sgn 


^( 1 - 0-1 


M- 


21 


M. 




i->-0 


7W- 


(48) 


In Eq. (l4^ I must satisfy I <C I but must not be so small as to render the numerator of the difference quotient equal 
to zero for too many lines in the ensemble. The systematic error introduced by this linearization of a derivative is 
controlled by 1. We note that if either or Ai- should be zero then the i9sgn term vanishes. However, this can 
only happen if either y_ or y_|_ is equal to zero. Although such worldlines exist, their number decreases with increasing 
Uppi, i.e., they form a set of measure zero in the worldline integral. 


D. N^umerical results for 7oo(xcp,t) aiid Tzz{x cp, t) 

Most sources of errors are already known. The statistical error is related to the number N of lines in our ensemble 
and the number of points per line Uppi determines a systematic error of the Monte Carlo integration. Furthermore, 
hs controls the systematic error of the derivatives (±y±)' and (±y±)". Since the worldline ensemble for closed loops 
is spherically symmetric in the continuum limit, the values of ((—?/_)) and ((?/+)) are identical and so are their error 
estimates. Therefore, we can immediately use the optimized values of these three parameters that we found for the 
single plate configuration. 

The linearization of 9 sgn in d^Ai introduces yet another systematic error that is determined by L In order to find 
an optimal value of I, we compute the expectation value of j{zcp/a) = dx sgn(x — 1) with x = AI+/AI_ as a function 
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a a 




a a 

FIG. 10. Worldline average of the derivative of the sign function in between the plates for different values of the linearization 
parameter 1. Fluctuations and errors are smallest for I = 0.075. 


of Zcpja. For that, we evaluate the difference quotient in Eq. (l48l) for different values of ^ <C 1. We show the results 
in Fig. [ini The function j{zcpla) is positive with large values around Zcp = 0 and decreases to zero near the plates. 
j{zcplo.) is only non-zero, if a worldline switches from intersecting one plate to intersecting the other plate first. This 
happens more often when the worldlines are in the middle than when they are close to one plate. For I — 0.1 and 
I — 0.075 the average of the derivative becomes a smooth function with small statistical errors. For smaller I the 
errors increase and the fluctuations around the average value become large. This is due to the fact that for too small 
I the numerator of the difference quotient tends to zero for too many worldlines in the average. We choose I = 0.075 
as an optimal value. The remaining errors are reduced by an increase of Uppi and N. 

We ascertain the overall systematic errors by determining how Too(afcp, t) and T^zixcp, t) converge to their analytic 
values for increasing Uppi. In Fig.[TT]we show numerical values for one exemplary component of the energy-momentum 
tensor for three different Zcp, using ensembles with N = 25 ■ 10^, Uppi G {10,..., 20} and the corresponding optimal 
value of hs and I = 0.075. We again fit our data to a function a+b/uppi such that a estimates the remaining systematic 
error for Uppi —>■ oo. We display those error estimates in Tab. IIIII The systematic errors vary with Zcp and so we take 
the largest error, 2.5%, as a general estimate. To arrive at results with minimal errors we increase the number of 
loops to = 5 • 10® and the number of points per loop to Uppi = 2^°. This value corresponds to / = 12 and e = 0.099 
according to Tab. ID We also choose I = 0.075. Fig. IT^ shows our optimized numerical data in comparison with the 
analytical result for the null energy condition in d = 2 and d = 3 dimensions. In both cases the null energy condition 
is violated at every point in between the plates with divergences near the plates. The differences between analytical 
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FIG. 11. Numerical values of Tzz\j^ in d = 2 dimensions at Zcp £ {0.15, 0.0 ,—0.35} as a function of Ifuppi. The statistical 
errors are represented by the error bars, and systematic errors are given by the difference between the fits to our data (solid 
lines) and the analytical values (dotted lines) in the limit Uppi oo. The plots for the other EMT components, as well as those 
for d = 2>, are qualitatively the same. 





to 



CO 



Zcp 

0.15 

0.0 

-0.35 

0.15 

0.0 

-0.35 

A«y4%] 

Too{xcp,t)\^ 

1.2 

0.7 

0.6 

1.8 

1.0 

0.8 

A,y4%] 

Tzz{Xcp, t)\j^ 

0.2 

1.2 

0.1 

0.6 

0.6 

0.4 

A,y4%] 

Tzz{Xcp, t) 1 Jj, 

1.5 

1.0 

1.4 

2.5 

0.5 

0.5 


TABLE III. Systematic errors /S.sys for the original Casimir setup, two parallel plates, with ensembles of = 25-10^ worldlines 
at Zap £ (0.15, 0.0,—0.35}. These values are derived as the difference between fits to our numerical data and the analytical 
values in the limit Uppi —>■ oo (cf. Fig. Hill . We use the largest error of 2.5 % as a general estimate. 


and numerical values are satisfactorily small with statistical errors smaller than 0.6 % and the general estimate of 
2.5% systematic error. 


V. CONCLUSIONS 


We have shown how the worldline formalism can be applied to local operators like the energy-momentum tensor of 
a scalar field that is coupled to a static background potential The general properties for computations of the 

EMT are similar to those of quantum actions or energies: the fluctuations are represented by spacetime trajectories 
of the fluctuating particles parametrized by an intrinsic propertime the total amount of which is a measure for the 
spatiotemporal extent of the fluctuation. Information about fluctuation-induced quantities is thus encoded in the 
geometric properties of the worldlines relative to the Casimir configuration, such as intersection properties of the 
worldlines with the plates. The essential difference for EMT computations is that in addition to closed trajectories 
also open trajectories have to be considered as a consequence of the point-splitting definition of the composite operator. 
While closed worldlines are sufficient for the energy component Tqo, the stress components also receive contributions 
from open worldlines, i.e., from the change of intersection properties upon opening a closed loop. 

In the present work, we have used a particular choice of the potential <j{x) to impose Dirichlet boundary conditions 
for massless field fluctuations on Casimir surfaces. Generalizations to softer properties of the plates and to a finite- 
mass scalar field are straightforward. Away from the Casimir surfaces, our construction guarantees that the EMT 
components are finite. Worldline numerics can then be used to calculate directly components of the energy-momentum 
tensor and check energy conditions. 

From a technical perspective, the use of a point-splitting method suggests employing common point lines in the 
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FIG. 12. Results for the null energy condition (NEC) for two parallel Dirichlet plates in 2 and 3 spatial dimensions. It is 
violated at every point in between the plates and diverges near them. The worldline numerical result (points with error bars) 
matches very well with the analytical result (solid lines). 


worldline average instead of center-of-mass loops for quantum actions. In contrast to the computation of effective 
actions, local operators do not only involve powers but also derivatives of cumulants of worldline variables. As any 
(propertime) lattice Monte Carlo calculation, our method comes with a statistical error controlled by the number N 
of configurations (loops or lines) in the ensemble, as well as with a systematic error from the discretization of the 
worldline propertime in terms of the number of points per loop Uppi. In addition, the approximation of derivatives with 
difference quotients leads to further systematic errors. We have demonstrated for several examples how to estimate 
and control theses errors straightforwardly. 

On the physics side, we have used our new method to re-calculate several results for the canonical energy-momentum 
tensor induced by a minimally coupled massless Dirichlet scalar. The well-known violation of the null energy condition 
is reproduced as well. Furthermore, we have shown how the conformal energy-momentum tensor can be constructed 
from its canonical counterpart within our approach. The conformal EMT also violates the NEC. Our numerical 
algorithm proved to be highly efficient, with results within at worst 3 % of the analytic values at computation costs 
of about I CPU-day on standard desktop PCs. We believe that our new method is now sufficiently mature to study 
more involved Casimir geometries in order to check the status of more sophisticated energy conditions such as, e.g., 
the ANEC along the lines of and beyond. 


Appendix A: Evaluation of lmG{x,x',k) 


In Sec. Ill Al we expressed the Green’s function G{x,x',k) in the worldline formalism. Due to the functional 
structure of the components of the energy-momentum tensor in Eq. (1151) . viz. T.. oc J dkf{k) Im J dsg{s) e*'*^ , the 
Wick rotations of s and k should not be performed independently but simultaneously. To do so, we pull the k integral 
inside the argument of the imaginary part. We use the Minkowskian versions of Eq. m and Eq. (HU) as well as 







21 


= m? + and thus find that only two different integrals (Ii and I 2 ) occur in Eq. (ITSl) : 


Ii = Im 


= Im 


i —ij dr 
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For both calculations we have performed the substitution s = —iT and k = —ikE in the integration variables and the 
rotation of the integral contours in the second line. We assume m > 0 and T > 0 for convergence. In the last step 
we have considered to be massless. The massless limit is, however, not a necessity and our calculations easily 

generalize to the case of finite m. It should be noted that the massless case is obtained as m —>■ 0“*', but can also be 
conceived on its own right by using Ei~ = '/W, which guarantees the correct analytic continuation under the double 
Wick rotation. 


Appendix B: Propertime integration of derivatives of the 0 function 


In Sec. Ill Cl we compute compact expressions of EMT components. This involves the interchange of several limits. 
Especially the propertime integration of derivatives of the step function deserve more details. These integrals can be 
solved by substituting q = y/TM. — 1 and integrating by parts. Then the integrals are of the following form: 


00 

A= [dq(-— ^ dqQ{q) + -- -d‘iQ{q) 

-1 


= lim 


r-^-l (g + 1)“ 1 




dq 
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_ (a-l)g 

((/+!)“ ^ (9 + 1 )“ 


5{q). 


(Bl) 


Here T is a function of Ai and its second derivatives and Q a function of M. and its first derivatives with respect to 
either 5z or Zcp- The exponent a is either d + 1 or d — 1. We identify dqQ{q) = S{q) since q = 0 is within the domain 
of integration. This makes the remaining integral trivial to evaluate. The boundary term vanishes because d{q) is 
zero everywhere except at g = 0 . 
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